AD  744508 


4*aMtaMMI 


*  V,  ' 


•'4  •*  -  . 

i. 


FOR 


t 


AEROSPACE  STUDIES  -  -  -  '  . 

•  •  ..  ;  ,  ■  •  *  -  ■*  ;r.*<  •>--** £  v. 


UNIVERSITY  OF  TORONTO 


tiflSk  *JR  -.72  si  2  7 & 


NKAR-FREE-MOLECULE  FLOW  CALCULATIONS- 
MONTE  CARLO  STUDY  OF  SKIMMER  INTERACTIONS 


N.A.  Derzko 


October,  1971 « 


'i . C  i  •'  •'  *  -'t  i 

v  • :  ;  .  ■  ! 

1  *  I..  - . — .-p-i 

UtllbLijJioLi  U  LlsiJ: 
•'  •;  A 


...JC'V 


t7  1  r, ,\v  1  jis:  i  \ 


R*produced  by  , 

NATIONAL  TECHNICAL  tfriA>  Technical  Note  No. 109 
INFORMATION  SERVICE 

U  S  Depoitmtnl  of  Comm*rc«  ; 

Spr.ngfield  VA  22151  ; 


1 


I 


I 


UKCLASSr  3D  J 

PAT 

’  Security  Classification 

DOCUMENT  dONTROL  DATA  -  R  &  D 

(Security  cttttiVet'loci  of  tlttm,  body  of  mbntrtct  mnd  indexing  mmoittion  mutt  bt  snterod  whtn  tht  ovtrmlf  rtpoet  it  cfdttHitd) 

1.  ORIGINATING  ACTIVITY  (Co* rotate  author) 

UNIVERSITY  OF  TORONTO 

INSTITUTE  FOR  AEROSPAfE  STUDIES 
■  TORONTO  ’5,  ONTARIO,  CANADA  ‘  i 

2a.  REPORT  SECURITY.  CLASSIFICATION 

.  UNCLASSIFIED 

S5.  <Jr6up 

3.  REPORT  TITLE  1  1 

NEAR-FREE-MOLECULE  FLOW  CALCULATIONS-MONTE  CARLO  STUDY  OF  SKIMMER  INTERACTIONS 

i 

4.  DESCRIPTIVE  NOTES  ( Typ a  ol  rape tt  and  incfu.tr.  da  tab)  t  : 

Scientific  Interim  i 

i  ' 

3.  AUT  NOFMSt  (F>:tt  nmmt,  miodlt  initimt,  Imtt  nmmt) 

i  *  » 

1  !  * 

N  A  DERZKO 

i  1 

|*  REPORT  DATE 

Pctober  1571 

7m.  TOTAL  NO.  OF  PAGES  176.  NO.  OF  REFS 

li  i  8 

•a.  CONTRACT  OR  GRANT  NO. 

AF0SR- 1481-68 

•  i  i 

b.  PROJECT  NO. 

9’^3-01  1  i  ; 

f-'6'.i02F 

681307 

«f. 

9a.  ORIGINATOR’S  REPORT  NUMBERIJ) 

UTIAS  T.N.  NO.  169 

i 

9b.  OTHER  REPORT  bUUSL&Sny  otharnumbaralhat  mar  ba  aaalgnad 

AroSB -1^72  *f2  78 

tO.  DISTRIBUTION  STATEMENT  i 

*  1.  ‘ 

Approved  for  public  release';  distribution  unlimited. 

1 1,  SUPPLEMENTARY  NOTES 

TECH,  OTHER  . 

i  ' 

12.  SPONSORING  MILITARY  ACTIVITY 

AF  Office  of  Scientific  Research  (NAM) 

1400  .Wilson  Boulevard 

Arlington,  Virginia  22209 

13.  ABSTRACT 


The  Monte  Carlo  simulation  'method  of  G.  A.  Bird  was  applied  :to  the 
.study  of  flow  through  a  skimmer  using  various  geometries,  speed 
ratios,  and  Knudsen  numbers.  The  computer  program  followed  a 
sample  of  3500  hard  sphere  molecules,  accumulating  distribution 
function  profiles  in  each',of  25  regions  of  the  flowfield,  tne  beam 
profile  at  tne  downstream  end  of  the  flowfield,  and  the  mass  flow 
through  the;  skimmer  orifice.  The  results  show  beam  intensities  of 
the.  order  of  .3  of’ the  no  skimmer  interaction  value  at  Knudsen 
numbers  of  01  and  demonstrate  clearly  tnat  most  of  the  interaction 
occurs  from '0  to  4  mean,  free  paths  upstream  of  the  orifice.  ' 


DD  ,Fr„1 473 


la 


UNCLASSIFIED 


,  Security  Classification 


A 


UNCLASSIFIED 


rant ty  CUiiliiotlon 


SKIMMER  INTERACTION. 
MONTE  CARLO  SIMULATION 
BOLTZMANN  EQUATION 


MOLECULAR  BEAM 


LINK  A 


HOLE  I  WT 


ih 


UNCLASSIFIED 

/  Security  CUtsificition 


NEAR-FREE-MOLECULE  FLOW  CALCULATIONS- 
MONTE  CARLO  STUDY  OF  SKIMMER  INTERACTIONS 

by 

N.  A.  Derzko 


I 

t 

i 

♦ 

i 

Submitted  March,  1971. 


i 


October,  1971. 


ITT  IAS  Technical  Note  No.  169 


ii 


aCKL’TOWLEDGEMENTS 

I  am  indebted  to  Dr.  G.  N.  Patterson  and  to  Dr.  J.  H. 
deLecuv  for  the  part  they  have  played  in  suggesting  and 
assisting  this  project,  and  to  J.  Kesten  for  programming 
assistance  in  its  early  stages.  This  work  was  supported 

by  the  U.  S.  Air  Force  Office  of  Scientific  Research  under 
Grant  No.  AF-AF0SR-68-1481A.' 


iii 


SUMMARY 

The  Monte  Carlo  simulation  method  of  G.  A.  Bird  was  applied  to  the 
study  of  flow  through  a  skimmer  using  various  geometries,  speed  ratios,  and 
Knudsen  numbers.  The  computer  program  followed  a  sample  of  3500  hard  sphere 
molecules,  accumulating  distribution  function  profiles  in  each  of  25  regions 
of  the  flowfield,  the  beam  profile  at  the  downstream  end  of  the  flowfield, 
and  the  mass  flow  through  the  skimmer  orifice. 

The  results  show  beam  intensities  of  the  order  of  .3  of  the  no  skimmer 
interaction  value  at  Knudsen  numbers  of  01  and  demonstrate  clearly  that  most 
of  the  interaction  occurs  from  0  to  4  mean  free  paths  upstream  of  the  orifice. 
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INTRODUCTION 


This  report  deals  with  part  of  the  work  we  have  been  doing  during  the  past 
three  summer a  on  the  development  and  testing  of  a  number  of  computer  program a 
designed  to  simulate  molecular  flow  problems  in  the  transition  regime.  The 
simulation  is  done  with  hard  sphere  model  using  the  Monte  Carlo 'method 
developed  by  0.  A.  Bird  (2). 

We  describe  here  its  application  to  the  study  of  the  flew  through  a  skimmer 
under  conditions  similar  to  those  involved  in  the  production  of  molecular  beams 
of  high  intensity  and  speed  ratio.  We  simulate  the  flow  which  occurs  in  a 
cylindrical  region  reaching  upstream  and  downstream  of  the  skimmer  but  not 
including  the  nozzle. 

Since  its  proposal  by  K&ntrovitz  and  Grey  in  1951 >  the  free  Jet-skimmer 
combination  has  proved  far  superior  to  the  older  methods  using  the  effusive 
flow  from  an  orifice  in  an  oven  source.  Nevertheless,  in  practice  the  intensity 
of  a  beam  produced  in  this  way  is  substantially  lower  than  that  one  would  obtain 
if  the  skimmer  just  removed  the  portion  of  the  free  Jet  around  the  beam  region  - 
that  is,  did  a  perfect  Job.  Considerable  experimental  work  has  gene  into  studying 
the  nature  and  causes  of  this  effect.  Recently  Govers,  LeRoy  and  Deckers  (6) 
have  studied  the  effects  of  background  scattering  on  skimmer  interaction;  Bossei, 
Huribut,  and  Sherman  (3;  have  done  experiments  to  study  the  magnitude  of  the 
skimmer  interaction  for  various  geometries  and  jets. 

The  effect  of  the  skimmer  interaction,  whatever  its  nature,  is  quite  clear 
If  one  plots  beam  intensity  in  particles  per  sec.  per  steradian  vs  r.czzle  skimmer 
separation,  the  graph  starts  at  a  low  value  fer  small  separation,  increases  to 
a  maximum  and  then  drops  rapidly  as  the  separation  exceeds  the  distance  to  the 
Mach  dxsc  (Fig.  4,  p.  99U ,  6 ) .  More  or  less  successful  attempts  have  been  made  to 
correlate  the  large  amount  of  experimental  data  available  to  enable  one  to  describe 
a  given  situation  with  a  fairly  small  number  of  parameters  With  skimmers  of 
similar  geometry,  the  ratio  of  mean  free  path  to  skimmer  orifice  diameter  appears 
to  be  significant  Roger  Campargue  [k)  has  discussed  a  number  of  dimensionless 
quantities  relevant  to  the  nozzle  beam  apparatus. 

The  theory  associated  with  a  nozzle  beam  divides  naturally  into  Three  areas: 

(a)  Description  of  the  free  Jet  -  that  i»  the  density,  velocity, 
and  speed  ratio  at  every  point  downstream  of  the  nozzle 

lb)  SKimmer  interaction. 

(c;  Development  of  the  beam  downstream  of  the  skimmer  (assuming 
that  any  molecules  involved  in  skimmer  interaction  leave 
the  beam) . 

A  satisfactory  theory  for  (a)  was  published  by  nshkanao  and  Sherman  in  Oft- 
(1).  A  number  of  authors  have  also  adequately  studied  be*,  eg.  Zapata  at  *i 
m  I960  ^3).  A  summary  appears  in  Hagena  and  Mar' on  (7).  Thus  far  (o'  has 
eluded  a  satisfactory  theoretical  discussion  it  the  purpose  this  vo.$  • * 
shed  some  light  on  this  dark  area  in  the  theory  of  noezie  beams. 


DESCRIPTION  OF  THE  PROGRAM 

? 

51.  METHOD  USED 

i 

The  Monte  Carlo  method  used  in  this  study  of  skimmer  interaction  is  basic¬ 
ally  the  method  of  G.  A.  Bird  which  has  been  described  by  him  in  a  number  of  , 
publications,  for  example  (2).  We  Include  a  short  description  here  for  sake  of 
completeness  and  to  enable  us  to  explain  the  changes  ve  have  Introduced  to 
accommodate  the  peculiarities  of  this  problem. 

• 

Bird's  method  is  a  Monte  Carlo  Algorithm  for  constructing  a  statistical 
sample  of  given  size  N  of  the  the  molecular  paths  occurring  in  sbme  finite 
regions  of  an  actual  flow.  The  construction  process  is  an  imitation  of  vha^, 
occurs  in  nature,  at  least  insofar  as  probability  distributions  are  concerned. 
The  particle  paths  are  functions  of  time,  so  that  at  any  instant  only  the  H  : 
coordinates  in  phase  space  (position,  velocity-space)  appear  In  computer 
memory.  A  time  counter  T  is  also  kept  in  memory. 

During  each  iterate  of  the  calculations  T  is  advanced  by  ,&  preset  quantity 
DTM  (Actually  a  number  DT  approximately  equal  to  DTM)  and  the  coordinates  in 
the  molecule  list  are  updated.  . 

§2.  INITIAL  CONDITIONS 

•  t 

The  initial  sample  is  from  a  distribution  based  on  what  we  would  have  if 
the  skimmer  absorbed  every  molecule  which  hit  it.  Actually,  the  initial 
distribution  is  of  secondary  importance  since  experience  has  shown  that  the 
molecule  list  approaches  the  correct  steady  state  distribution  in  relatively 
few  iterations. 

As  soon  as  T  reaches  a  preset  value  TEQ  (time  at  which  steady  flow  is 
assumed  to  have  been  achieved)  data  collection  for  subsequent  output  i3  begun 
and  continues  until  T  reaches  TEND.  Following  output  the  program  automatically 
stops  or  does  the  next  case.  i 

% 

53.  ONE  ITERATION 

Each  iteration  has  two  components: 

(a)  Changes  to  molecule  list  due  to  collision.  ; 

; 

(b)  Changes  to  molecule  list  due  to  convection. 

The  computation,  of  (a)  requires  that  the  flow  region  be  divided  into 
cells  (Fig.  2)  in  such  a  way  that  the  distribution  function  remains  fairly 
constant  over  each  cell.  At  any  moment  the  molecules  residing  in  a  given 
cell  constitute  a  sample  from  the  distribution  for  that  cell.  If  these 
molecules  are  allowed  to  collide  among  themselves  then  these  collisions  are 
representatives  of  the  collisions  actually  taking  place  in  the  given  ce.V'  The 
theory  backing  this  procedure  is  stronger  than  the  above  suggests.  We  refer  to 
Bird's  original  papers  on  the  subject. 

* 

The  details  of  the  collision  process  are  as  follows:  A  pair  of  molecules 
is  chosen  in  such  a  way  that  the  probability  of  being  chosen  is  proportional  » 
to  the  relative  velocity.  If  and  ?h  are  the  velocitiea  of  the  pair  choaen* 
they  are  replaced  according  to  the  hard  sphere  collision  law  with  v%,  where 
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where  n  is  ’chosen  from  the  uniform  distribution  on  the  unit  sphere. 

•  1  *  , 

The  position  poordinates  rema'in  unchanged  and  a  time  parameter  D T  is 
advanced  by  .  i 

1  1  I 

I 

,  DDT  «  TCONS!  « 

“2|Vval 

,  The  number  TCONS  is -fix'ed  fo*  each  tell  and  is  explained  in  £6.  Collisions 
are  computed  for!  a  given  cell  until  ,  ! 

.  ,  •  i' 

s'  DT  +  DDT  >  DTM.  1  •  : 

'  -  ■  ' 

The  last  collision  is  kept  with  probability 

'  -  i 

DTM1  -  DT  • 

'  . .  DDT 


and  ignored  otherwise.  This  last  refinement  is' essential  here  because  some  cells 
have  isu-'h  low  collision  frequencies  that  no  collisions  at  all  should  be  computed 
during  -some  intervals  DTM.i 

*  i  ' 

,  The  computation  of  (b)  is  basically  simple,  the  position  coordinate^  are 
translated  on  the  basis  of  the  time  interval  DTM  and  velocities.  The  computation 
involves  some  wifirk  because 'we  use  cylindrical  coordinates  in  the  program  in  which 
the  translation  equations  are  nonlinear We  have  had  to  introduce  two  additional 
complications  which  do  not  appear  necessary  at’,  fir's t  sight  .—  referencing  the 
molecule  list  §9'  -  and  stock  Splitting  §5.  At  the  end  of  the  collision  convection 
■cycle  T  is  'incremented  appropriately  and  the  whole  process  is  repeated. 

•  .  i 

■§U.  i  GEOMETRY  , 

,  i 

Sines  the  skimmer  problem  has  cylindrical  symmetry  about  -.the  common  axis  of 
the  free  jet  and  the  skimmer,  twp  coordinates  suffice  toi specify  position  -  the 
axial  ahd  radial..  They  are  denoted  by  (x,y),  in  the  program  (rather  than  \o,r) 
as  Is  customary  because  essentially  the  same  program  with  a  few  changed  subroutines 
car^  be  used  ih  Cartesian  coordinates ) ,  The  velocity  coordinates  (u,v,w)  are 
with  respect  to  the  Cartesian  ftfame  I,J,k  wher,e  I  is  in  the  axial  direction,  3 
is  in  the  radical  direction  and  k  is  in’  the  direction  of  increasing  0.  Tne  actual 
region  of  the  flow  under  study  ,1s  the  ^olid  of  revolution  ^generated  by  the  shape 
in  Fig.  1.  |  i  ' 

The  flow  region  is  subdivided,  into  cells  as  requited  in  the  collision  cal¬ 
culation.  The  25  cells  and  the  parameters  necessary  to  specify  the  geometry 
are  marked  in 'Fig.  1.  These  parameters  are  kept  in  common  area  (GEOM)  (see 
Appendix  1).  Some  of  them  must  be  specified  initially  (see  Block  Data  subprogram' 
and  the  rest  are  competed  'in  subroutine  SETUP.  ,  , 

STOCK  3PlllTTING  1  ’  ’  , 

1  !  ‘ 

'  If  the  coordinates  in  computer  memory  simply  constituted  a  sample  of  the 
molecular  coordinates  in  an  actual  flow,  we  should  find  'that  their  density  in 
the  x.,y-plane  ’  increased  linearly  with  x.  This  would  be  unsat  1'sfactory  because 

then  most  of  the  computation  would  be  i  done  for  the  region  far  from  thb  axis  vtiue 

1  '  ) 


-li¬ 


the  statistical  scatter  for  the  beam  region  containing  only  a  small  number  of 
molecules  would,  render  these  results  almost  useless. 

The  solution  we  have  adopted  to  this  difficulty  involves  subdividing  the 
flow  region  into  a  number  of  subregions,  using  the  streamlines  of  the  flow  to 
construct  the  boundary  surfaces  whenever  possible.  The  subregions  are  then 
assigned  weights  in  such  a  way  that  a  simulator  molecule  in  a  subregion  far 
from  the  axis  ’’represents"  several  times  as  many  "real"  molecules  as  one  in 
an  axial  region.  For  programming  reasons  we  have  found  it  convenient  to  make 
the  weights  inversely  proportional  to  the  representation  multiples,  so  that 
the  subregions  farthest  from  the  axis  have  the  lowest  weights. 

Such  an  arrangement  calls  for  a  special  strategy  when  a  molecule  happens 
to  move  between  subregions  of  different  weights.  The  strategy  is  this:  When 
a  molecule  moves  from  a  weight  w^  region  to  a  weight  W  region,  it  is  replaced 
with 


INTR(*£/'wi) 

identical  molecules.  The  random  variable  INTO  is  defined  as  follows: 

[a]  with  probability  l-o+[a] 

INTO  (a)  = 

[a]  +1  with  probability  a- [a] 

where  [a]  denotes  the  greatest  integer  less  than  or  equal  to  a  as  is  customary. 
We  have  called  this  process  "stock  splitting". 

56.  TCONS  CALCULATION 

The  original  formula  of  G.  A.  Bird  for  the  time  advance  per  collision  is 
DDT  =  2fij?/(&|jij?v)  where  0^,  n^,  a^  are  respectively  the  volume,  total  number  of 
simulator  molecules,  and  cross-section  per  simulator  molecule  for  cell  jc,and  v 
is  the  relative  velocity  of  the  collision  pair.  When  the  cells  are  weighted 
then  a^  must  be  inversely  proportional  to  the  weight.  That  is  if  the  weights 
are  wK,  then 


WK  =  COnSt  =  ** 

The  additional  requirements  on  in  an  actual  calculation  are 

(1)  The  cross-section  density  in  a  specified  region  of  the  flow 
is  given  (see  §10). 

(2)  The  total  number  of  molecules  must  be  in  the  vicinity  of 
some  Nq  smaller  than  the  capacity  of  the  molecule  list. 


Let  us  say  N0  is  the  number  of  molecules  obtained  in  the  flowfield  if  the 

skimmer  caused  no  disturbance,  i.e.  absorbed  everything  that  hit  it.  Let  AK 

be  the  total  cross-section  ( (cross' -section/molecule)  x  number  of  molecules)  in 

cell  t.  A.^  can  be  computed  from  the  free  jet  structure  and  the  cross  section 

density  A  at  the  skimmer  orifice  (subroutine  VPOL).  Then 
o 


N. 

< 


-  fjc  -  Vfi 

*K  '  X" 


T 


and 


so  that 


and 
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^  A*wk: 


TC0N8(k)  =  2^/a^.  =  2J1.  w  N0/IA  w 


DDT 


TC0N5U; 

2 

n  v 
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§7.  THE  INCOMING  DISTRIBUTION 

First  the  total  influx  rate  CU  to  the  flow  region  from  the  jet  is  computed 
from  the  jet  cross-section,  density  Ae  at  the1 entry  point,  its  velocity  v  (taken 
=  1  in  our  units)  and  the  weights  of  the  cells  receiving  the  molecules.  The 
total  entering  cell  k  (assuming  all  k  borders  on  the  entry  plane)  is 

A  v 
e 

-  x  area  of  entry  surface  element . 

< 

CU  is  the  sum  of  all  such  quantities. 

The  number  entering  the  fZ  -w  region  during  time  DT  it  CU  x  DT.  Each  entering 
molecule  has  its  position  coordinates  chosen  at  random  m  accordance  with  the 
Ashkenas-Sherman  formula  [l]  for  variation  of  jet  density  with  distance  from  the 
axis 


p(R,o) 


and  the  weight  of  the  ceil  about  to  receive  the  molecule.  The  details  cl  this 
calculation  appear  in  subroutine  INPUT  (-Appendix  1), 


§8.  BOUNDARf  CONDITIONS 


The  boundaries  in  question  are  the  skimmer  walls  and  the  artificial  boundaries 
introduced  to  enclose  the  flow  field.  A  molecule  colliding  with  a  skimmer  has  a 
choice  of  being  reflected  or  reemitted  from  a  cosine  distribution  at  a  specified 
temperature  CR.  The  probability  of  reemission  is  a  specified  accommodation 
coefficient  ACCOM. 

Tne  artificial  boundaries  are: 

(1)  The  entry  plane, 

(2)  The  upstream  outer  boundary, 

(3)  The  downstream  boundary. 

Molecules  which  collide  with  (1)  and  (3)  are  deleted.  If  a  molectue  collides  with 
(2)  its  velocity  vector  is  reflected  and  augmented  by  the  vector  2pn  where  r.  is 
the  outward  normal  to  the  wall  the  p  is  a  specified  constant.  If  the  resulting 
velocity  vector  points  out  of  the  region  the  molecule  is  deleted;  otherwise  it 

\ 
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is  allowed  to  reenter.  When  p  is  roughly  equal  to  the  normal  component  of 
the  average  velocity  of  the  molecules  crossing  the  wall,  its  net  effect  is 
to  simulate  no  discontinuity  in  the  flow  due  to  the  boundary. 

§9.  REFERENCING  THE  MOLECULE  LIST 

We  discuss  here  the  actual  arrangement  of  the  molecule  coordinates 
C^  d  *  1,2,... in  computer  memory  consistent  with  a  fast  method  for  referring 
to  the  coordinates  corresponding  to  a  given  cell. 

The  list  is  arranged  in  such  a  way  that 

C1  *  *  *  CN1 


are  in  cell  1, 

CN1+1  '•*  CN2 

are  in  cell  2,  etc.  The  numbers  NI  are  stored  in  an  auxiliary  li3t  called  MAP. 
The  basic  operations  which  must  be  performed  on  this  list  are 

(1)  Removal  of  a  molecule,  and 

(2)  Addition  of  a  molecule. 

(1)  To  remove  C*  where  Nj.,  <  <  <_  N, ,  store  C  in  place  of 
place  of  C^ ,  etc.  until  the  end  of  the^list  is  reached. 

(2) .  To  insert  coordinates  C  -into  cell  i  we  make  room  for  it  by  storing 

the  first  molecule  of  each  cell  in-' the  location  following  the  last  molecule, 
starting  from- the  end  and  going  up  to  cell  i  +  1.  This  leaves  location  N.  +  1 
for  C.  1 


c*’  cyut  iD 


Of  course,  in  each  case  MAP  is  updated  to  reflect  the  new  configuration. 

It  is  clear  that  a  simple  transfer  of  molecule  from  one  cell  to  another 
does  not  require  the  operations  to 'go  to  the  end  of  the  coordinate  list.  A 
time  saving  procedure  to  accomplish  this  task  is  incorporated  into  the  program. 

§10.  SIMULATING  AN  ACTUAL  EXPERIMENT 

The  program  takes  the  limiting  mean  centreline  jet  velocity  to  be  1.  All 
other  velocities  are  given  in  terms  of  it.  The  simulator  geometry  must  be 
similar  (in  the  geometric  sense)  to  the  real  geometry,  but  the  units  we 
arbitrary.  Once  the  velocity  and  the  distance  unit  are  fixed,  the  time  unit 
is  automatically  determined.  Since  we  are  mostly  concerned  with  the  steady 
state  we  do  not  even  need  to  compare  computed  times  with  real  times  except 
in  the  initial  approach  to  the  steady  state.  '  , 

,The  remaining  parameter  that  .->ust  be  assigned  is  the  cross  section  density 
A  at  some  reference  point  .  the  flow  which  we  have  chosen  to  be  the  skimmer 
orifice  assuming  perfect  skimming.  A  has  units  of  reciprocal  length,  in  fact 
it  is  Just  a  multiple  of  1/mean  free  path.  Thus  Ad,  where  d  is  the  skimmer 
orifice  diameter  is  dimensionless,  and  we  give  it  the  same  value  in  the 
simulator  as  it 'has  in  the  experiment. 

Though  this  prescription  is  simple,  it  has  its  problems,  not  the  least 
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being  that  the  hard  sphere  cross  section  of  real  molecules  depends  on  their 
relative  velocity.  The  best  we  can  do  is  to  use, some  sort  of  average  relative 
velocity.  The  realtive  velocity  of  the  beam  and  the  molecules  coming  off  the 
skimmer  wall  is  a  good  candidate  because  ;his  is  when  most  of  the  collisions 
take  place. 

511.  OUTPUT  FORMAT 

The  output  from  each  run  contains  the  following  information  (see  Appendix 

2). 

(1)  A  printout  of  DT  and  MAP  foi  each  interval  DT'M  for  checking 
purposes . 

1 2 )  A  flow  density  profile  at  the  skimmer  orifice  as  well  as  a 
total  particie  flow  rate. 

(3)  A  flow  density  profile  at  the  far  downstream  end  of  the 
flowfield  as  well  as  the  total  particle  flow  rate  through 
the  geometrical  beam  area.  This  last  figure  is  useful  m 
estimating  the  centreline  beam  intensity 

(4)  Histograms  for  the  integrated  distribution  function  profiles 
/ fdvdw,  / fdwdu,  / fdudv  for  each  cell. 

The  distribution  function  is  broken  down  into  the  distribution  function 
of  the  molecules  which  have  not  undergone  collision  in  the  xlowfieid  and  the 
remainder.  The  mean  velocity,  mean  square  velocity,  and  mean  number  of  moiecules 
are  also  given  along  with  a  theoretical  comparison  number  for  no  skimraev  inter¬ 
action.  The  histograms  are  obtained  by  time  averaging  regular  observations  of 
the  flowfield  for  the  duration  of  the  computation  after  equilibrium  is  reached 

RESULTS 

A  simulator  program  such  as  tms  one  should  be  viewed  rather  i_k.-  an^ 
experimental  facility  suitable  for  investigating  a  range  of  problems  connected  .. 
with  the  sximmer.  We  describe  here  tne  results  of  a  study  of  SKimmec  interaction. 

A  preliminary  study  was  done  to  checx  a  modex  proposed  by  French  ..  3)  suggesting 
a  significant  interaction  downstream  of  the  skimmer  oriir.e.  This  was  done  in  a 
program  which  preceeded  the  present  ono  by  making  the  downstream  oKimmer  waii 
absorbing.  No  downstream  effect  was  observed.  A  check  of  the  distribution  functions 
m  the  present  study  supports  this  con-iusion. 

In  fact  this  study  shows  that  the  skimmer  interaction  is  caused  by  a  screen 
of  molecuies  originating  from  the  wall  and  from  collisions  extending  approximately 
one  mean  free  path  upstream  oi~the  skimmer  orifice  Tne  supporting  materia*  for 
this  conclusion  is  contained  in  the  following  grapns  ana  tables 

Table  i  shows  the  cases  computed  m  this  study. 

Fig.  i  defines  the  geometric  variables  in  the  program  ana  Fig's  2  -  give 
the  specific  geometries  used  in  the  calculations. 

*  '•  >. 

Table  2  shows  condition  near  the  skimmer  wall  fox  an  assortment  of  cases. 

*t  is  most  significant  that  for  ail  but  the  sharpest  skimmer  .angles  50-', 
the  flow  is  subsonic.  We  suspect  there  is  always  a  subsonic  region  ci-se  to 


"vrrrw 
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,  « 

j  ,  > 

the  skimmer  wall  which  is  very  thin  for  small  angles  but  grows  rapidly  as  ■  ■ 

the  angle  increases.  Also  the  surface'  interaction -with  the  .skimmer  plays 
a  significant  role,  since  case  23  (ACCOM  «  0)  gives  a  high  Mach! number  even 
for  a  blunt  skimmer.  •  ; 

*  ,i 

‘  ■  -r  I 

Fig.  5  shows  theimass  flow  rate  through  the  skimmer  opifice  vs  A  in  . 
terms  of  the  theoretical  no  interaction  value.  ,  The  mass  flow  is  less  than  • 
one  in  all  cases  computed  but  shows  a  tendency  to  go  to  1  in1  the  c8ld  t 

skimmer  and  no  lip  cases.  1 

i  (  , 

Fig.1  6  shows  beam  ixitensity  vs  A.  The  importance  of  the  skimmer  lip  '< 
sharpness  is  here  clearly  indicated.  The* skimmer  angle  seems' to  have  less 
effect  than  the  skimmfer  temperature.  Though  the  general  shape  of  the  exper¬ 
imental  beam  intensity  vs  Knudsen' number  curves  is  evident,  our  results  do' 
not  show  decreases,  in  intensity  as  strong  51s  those  in  [3]  for  instance.  «  1 

*.  |  S  * 

Figs.  7,  8  show  mass  flow  and  beam  intensity  vs  various  parameters.  It 
is  most  interesting  that  an  increase  in  PRES  which  leads  to  a  decrease  in  the, 
number  of  molecules  entering  the  flowfield  through  the  outer  upstream  boundary 
leads  to  a  decreased  mass  flow  but  no  significant  change  in  beam  intensity. 
Furthermore , ( a  low  accommodation  coefficient  gives  high  values  in  both  graphs. 

Figs.  9,  10  were  obtained  by  smoothing  the  histograms  produced  by  the 
program.  It  is  cle'ar  that  a  Knudsen  number  of  1/2  (A,  =:U)  lejads  to  more  or 
less  complete  destruction!  of  the  beajp  by  skimmer  interaction. 

»  . 

Figs.  11,  12  show  clearly  a  sharp  increase  of  molecules  having  a  bread  1 
distribution  function  and  infiltrating. the  beam  region  between  Knudsen  numbers 
1/2  and  1.  ,  1  '  ,  ,  . 

1  !  1  .  1 

Finally  .a  note  on  ihe  standard  'deviation  of  .he  results.  The  beam  intensity 
and  mass  flow  curves  have  a  relative  standard  deviation  of  between  3%  and  3%’ 

It  follows  from  the  theory  of  Monte  Carlo  simulators  that  the  probability  distri¬ 
bution  of  each  value  obtained  is  normal,..  Accuracy  is  expensivei  in  these  calcul¬ 
ations.  A  folxfold  increase  in  computing  time  would  cut  the.  standard-  deviation 
in 'half.  1 


'  1 
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'TABLE  1 

CASES  COMPUTED  • 


CHARACTERISING  PARAMETERS  '  _ COMMENTS 

A  SIA  .  CR  :  PRES  ACCOM  .THLIP  '  GEOM  ■  ' _ _ 


1 

;  6.5 

O.’l 

1.' 

1. 

0.5 

0.1v 

'3 

Sharp  skimmer 

2 

.1. 

0.1 

1. 

1-  : 

0:;5 

0.1  k 

3 

3  ' 

2. 

0,1 

1. 

1. 

Ov.5 

Oil 

!  3  ■ 

4 

4. 

091 

1. 

1. 

0'.'5  , 

O-fl 

3 

1 

5 

2. 

0.1 

1'.  • 

2. 

0.5 

0.1 

1 

Effect  of  PRES 

6 

2. 

0.1  1 

1. 

1.5 

0-5  , 

0.1 

1 

7 

2. 

.  0.1 

lrf 

0;75 

0.:5‘ 

0.1 

1 

8 

2.  . 

0.1 

1. 

0.5 

0.5 

0.1  1 

i  , 

* 

9 

2. 

0.0546 

1. 

1.- 

0.5 

:  0.1  . 

1 

High  mach  number 

10 

1. 

0.1 

1. 

2.  , 

015 

0.-1 

1 

Absorbing  walls 

11 

.  0.5 

0.1  ' 

0.5 , 

1.- 

0.5 

0..1 

1 

dold  skimmer  series 

12 

'  1.  ' 

0.1 

0.5 

1. 

0.5 

0.1 

■  ‘i 

13 

2. 

0.1 

0:5 

1. 

0.5 

O-.l 

1 

l4 

!  k. 

q.l 

0.5 

1. 

0.5 

0.1 

1. 

If 

0.5 

0.1 

>1. 

1. 

0.5  ! 

0.1 

1 

Standard  series 

16  - 

1.0 

0.1 

1. 

i. 

0.5 

0.1 

1  ! 

blunt  skimmer 

17 

2.0 

0.1 

'  1. 

1. 

0.5 

o.i 

1 

l8 

4;0 

•  0,1 1 

1. 

1. 

0.5 

0.1 

1 

'19 

0.5 

0.1 

,1.  ! 

1. 

0,5 

0.0 

1 

Perfect  edge . 

20 

1.0 

0.1 

1. 

1. 

0.5 

0.0 

1 

21 

!2.0 

0.1 

1. 

,  i-  ’ 

0.5 

0.0  • 

1 

i 

22 

k.O 

0.1  * 

1, 

1, 

0,5 

0.0 

i  ■ 

• 

23 

2. 

0.1  * 

1. 

1. 

0. 

0.1 

.  1  • 

Effects  of  surface 

2k 

‘  2.  1 

1  0.1 

1. 

1. 

1. 

o.i' 

1 

interaction 

25 

0.5 

0.1 

1. 

1. 

0.5 

0.1 

2 

Intermediate  angle 

26 

1. 

.0.1 

1. 

1. 

0.5 

0.1 

,  2 

27  ■ 

2. 

0.1 

1.' 

1. 

,  0.5 

0.1 

1  2 

28  • 

4. 

0.1 

: 

•  l. 

i. 

0.5 

• 

’  0.1 

1 

2 

i 


Notes 

A  vp 

Knuds en 

numbdr  Kn 

A 

0s.  5 

12  4  ! 

1 

Kn  : 

4 

2  1  »  0.5 

SIA,  CR,  PRES,  ACCOM,  THLIP  are  explained  in  Appendix  -• 
GECA  1,  2,  3:'are  shown  in  Figs.  1,  2,  3  respectively. 
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TABLE  2 

CONDITIONS  IN  FRONT  OF  SKIMMER 


CELL  14  CELL  15 


c*ss,. 

Mean  Speed 

Mach  # 

Mean  Speed 

Mach  # 

Inclination 

15 

0.55' 

0.64 

0.51 

0.55 

0.83 

18 

0.52 

0.79 

0.34 

0.44' 

0.69 

5 

0.62 

0.88 

0.48 

0.57 

0.52 

8 

0.33 

0.44 

0.28 

0.34 

0.73 

1 

0.80 

i.4 

0.72 

1.2 

0.28 

4 

0.86 

2.1 

0.71 

1. 

0.32 

23 

0.87 

2.2 

0.75 

2.3 

0.62 

25 

0.68 

1.0 

0.60 

0.82 

0.43 

28 

0.71 

1.3 

0.53 

0.73 

0.47 
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FIG.  1  DEFINITION  OF  GEOMETRICAL  P 


PIG 


FLOWFIELD  FOR  GEOM  =  3,  SKIMMER  ANGLES  50 


Normalized  Total  Flux 


Normalized  Beam.  Intensity 


FIG.  T  MAS3  FLOW  THROUGH  SKIMMER  VS  OTHER  PARAMETERS 


PRES  'S«r>9*  • 


FIG.  8  BEAM  INTENSITY  VS  OTHER  PARAMETERS 


PRES  Series  • 


Accom  *  I 


FIG.  11  DISTRIBUTION  FUNCTION  IN  CELL  4'  (GE0M*2) 
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COMMENTS 


A.  The  Variables 


Common  Area 


Al-1 


/RCRD/ 


7geom/ 

/NPRFL/ 


/CONS/ 


/BEGIN/, 


Store^e  area  for  beam  and  orifice  profile  histograms. 

BPRFL  and  HPRFL  resp.  DBPR  and  DHPR  are  the  widths  of  the 
histogram  celis  in  each  case.  ,, 

"  ,  ‘  •  . 

Geometrical  definition  of  the  flow  field  -  see  Fig. 

For  internal 'use  in  the  program  to  save  computer  time  in 
the  introduction  of  molecules  through  the  input  plane. 

1  CU  is  the  total  number  of  molecules  entering  the  flow  field 
through  the  input  plane  per  unit  time,  SIA  is  the  speed 
ratio,,  in  the'  jet,  CR  is  a  speed  defining  the  temperature 
.of  the  skimmer  walls,  ACCOM  is  the  accommodation  coefficient, 
PRES  is  relatecl  to  the  pressure  on  ithe  upstream  outer 
boundary, '  CliWT' is  the  list  containing  the  cell  weights,, 
TCtoNS  i£  the. list  .of  time  advance  per  collision  constants 
for  the  cells.  The.  remaining  variables  are  for  internal 
use  of  the  program  to  save  recomputing  frequently  used 
numbers.  They  are  initially  computed  in  SETUP. 

i 

,BE  and  §G  are  initial  values  fdr  the  sequencing  variables 
of  tne  random  number  generator^.  A  13  the  cross  section 
density  at  the  ’orifice,  and  NST  is  .starting  number  01 
simulator  molecules  m  the  flbw'  field. 


/TIME/ 


/STORE/ 


/HIST/ 


T  -  Accumulated  time  variable  for  the  program. 

DT  -  Weighted  time  per  iteration  (DT<^DTM). 

DTM  -  Preset  time  allowance  per  iteration. 

TEQ  -  Value  of  T  for  which  equilibrium  is  assumed  to  have 
been  achieved. 

TEND  -  Value .of  T  At  which  computation  ceases  and  output 
1  is  begun. 

Tne  remaining  variables  were  ndt  actually  used  in  this 
program J  -  !  ' 

NCL  is  the  number  of  celis,  MCP  =  5000,,  and  MAP  is  for 
purposes  of  referencing  the  list  of  molecules.  The 
remaining  variables  in  this  oomnpn  area  constitute  the 
molecular  list.  1  ! 

Storage  area  for  the  histograms  FU  etc;  CFU  etc;  the 
mean  velocity  vectors  for  each  c'ell  -  UB  etc.  and  tn= 
mean  square  velocity  vector  for  each  ceil  UUB  - 
CORG  counts  (the  number  of  molecules  in  each  ceil  which  have 
undergone  more  than  one  collision. 

Contains  control  parameters  for  the  data  Collection  operation. 
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/RAND/  Contains  the  sequencing  parameters  for  the  random  number 
generators . 

/VARS/  Temporary  storage  area  for  one  molecule  to  save  address 
computation  time  while  the  molecule  is  being  processed. 


B.  Subroutines _ 

MAIN  This  program  is  used  to  set  up  the  cases.  Each  call  to  PROG  constitutes 

the  computation  of  one  case. 

BLOCK  DATA  Used  to  assign  values  to  variables  which  are  common  to  the  computation 
of  many  cases. 

PROG  Control  program  for  the  computation  of  one  case. 

SETUP  Computes  variables  not  assigned  in  MAIN  or  BLOCK  DATA  and  performs  all 

the  initialization  necessary  to  begin  computation  of  a  case. 

VPOL  Computes  cell  volumes  and  cell  masses  for  the  undisturbed  free  Jet. 

INPR  Sets  up  an  array  used  to  obtain  the  input  profile  by  interpolation. 

LOC  Function  which  computes  the  cell  containing  the  coordinates  (x,y).  LOC 

works  by  computing  a  word  IX  which  has  0  or  1  in  the  k-th  bit  depending 
on  which  side  of  the  k-th  boundary  (x,y)  is  found.  Then  IX  is  compared 
with  a  standard  ordered  list  (MASK,  1DFR) . 

COLI  Collision  increment  for  cell  L. 

SQM  Estimates  the  maximum  relative  velocity  in  a  given  cell. 

TRBKP  "Translate  and  Book-keep".  This  subroutine  essentially  controls 

computation  of  the  convection  increment. 

INVRS  Computes  the  cell  number  L  of  the  M-Th  molecule  in  the  molecule  list. 

UP,  DOWN  Used  in  inserting,  deleting  or  moving  molecules  in  the  molecule 
list. 

Transfer  of  coordinates. 

Checks  a  given  molecule  for  collision  with  boundaries  and  it  moves  it 
a  distance  appropriate  to  the  time  DT. 

Changes  coordinates  of  a  given  molecule  appropriate  to  the  time  interval 
DT. 

LCHK  Checks  whether  collision  occurs  with  a  given  boundary. 

EMIT  Emits  a  fully  accommodated  molecule  in  the  direction  C,S. 

REFL  Reflects  a  molecule  in  the  direction  C,S. 

Introduces  CU  and  DT  molecules  through  the  input  plane  on  the  basis 
of  the  input  profile  defined  by  FNPR,  DFNPR. 


TRS,  TR 
FRP 

TRANS 


INPUT 
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C3SFL 

TABU 

FIT 

SIFT 

INTR 

RANDG 

RANDL 

RANDR 

UMAX 

ISGN 

DUMP. 

BAR 

HIST1 

HIST3 


I 

( 

# 

Cubic  Spline  interpolation  subroutine. 

Controls  the  accumulation  of  histograms  and  moments  of  the  distribution 
function. 

Inserts  molecule  at  cell  L. 

Auxiliary  subroutine  to  accumulate  histograms. 

See  S5. 

Standard  normal  random  number  generator. 

Uniform  on  [0,1]  random,  number  generator. 

Uniform  on  [-1,1]. 

Finds  the  maximum  element  in  an  array. 

ISGN(X)  ■  1  if  x  ^  0) 

0  otherwise 

Controls  the  output  function  of  the  program. 

Sets  up  one  histogram  bar  for  printing. 

Sets  up  one  histogram  for  printing. 

Sets  up  and  prints  the  3x3  histogram  describing  the  distribution 
function  in  cell  L. 
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APPENDIX  t 


PROGRAM  LISTING 


//FQNT.SYStN  Oft  •  — - 

71mo3SB  /ufAn/'nnkw.oMtHr.upBf n$o» 1 101 

common  /cfom/  pf l.nf *N*,m  an*. 1i*t  IP»%LQ«Sllt  *CM*xCB?«  vc*^**C"**  • 

x  *i  I  I.Si!?.xCAn,xCM./eM*vCAl.vCAJ.vcM.YCA4.tt*S.vC“N.vC*7. 
x  uni.sto/.MOx.sifH.sios.  siu»ui*.  raoin.oisnu  wuv»i, 

_ t  juxixiiiQx-ihiix»±-xnL  i  2  a  uaotn  i  zu - 

COMMOH  /NPMFI./  /NPR|  10)  »FNNH  IP)*  W"»N<  101 

common  /C on*/  cn.*1A, *IR,C»»CT, ACCOM, hRF S.CIMTI/m  .TCihsMM 

COMMON  /APr.tN/ftR,Mr.,A.N*T 
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VCRA»/.S*niAM* 

OI*N|.nt*N***CM 
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YV  f( XMJ  .YC.P7*  |  imviiiH) 

YVTI7?).YrA7*U  I74XV1I7T1 
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vm  M4|.vrfH.txv?,YVT,4,l?,i*,4*9*o.r»is"S.vwnii4)i 
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X  ♦  T.|4| sr.NISm*X«VCA?-V>  ♦  TK4|<,;n(SU?***vCA7“Y> 

X  ♦  T  l  •  I  SON  I XCA4-X  I  ♦  TM4dr.NI  If  AT-x  I  ♦  TNr  I  sr.N  |  IC  A?-X  } 

x  ♦  Tn4|sr.Mixi  *  TP4|sr.Ni«.xrAi i 
_ iMLJ.-J.Uj2S _ 

l|XH*| IX.4M0.1MASXI I ) 

I  I  -  t 

IF  lIXM.lnfAIIM  1*7*1 

1  C04TI4UF 
I0f*74 

-^1«M4 _ 

2  inc»ii 

AFTIM4 


suAnniiTiwf  cm  my 

U|«PFA  1970 


X  M(*000),iniSOtK>l 

COMMON  /fONS/  Cn.S|Ad|A,CA,CT*ArcnM,PAFS*CtWTI741,TCOM^(7M 

COMMON  /RANO/A.r. 

COMMON  /T I Mf/  T.OT »njMv  TF0*nfSM, TFNO.OTRP. TAFP 

•‘I  AM  FI  IST/rHFX/  XA.KA.MA.MA.t  f  FRV*||A.VA.MA.|IA*VA*NA*MAP 

_ _ 


IF  ll-ll  1.1.7 

«A«0 

Ml  TO  3 

KA«MAPIt-t  1 

A4-MAPI  it 
F*4A-Fl.OAT  IJA.kA) 


AA*4A«l 
IF  10.10.4 

frT.TfnNSIl»/F«AA/FxA4 

VNAX.^ORTI  S0«|  * A,«R,n»*SOM|XA,X*. y J ♦ V0*M  X A.X A , M M 

I F  |  VM At |  10.10.4 

MA.^A.  |Ff*(F*4A»rf  JNm  litn 


MH-KA*lF!t|F<MA4RANn(.(M  M 
IF  INN.MA)  4.4.4 
|IA-ll(  4A| 

VA«V(MA| 

R1,Y| MAI 
|iA*ii|«*A} 


V«4\ l -*| 

***N(M|  * 

U*A-|l4«MA 

V*AaV«-V A 

W4A.hH.WA 

«Fj  .4Q4T  m»  A»n*  A*  V4A4VAA.UR  A4WAA 


JFCVRFl  ,|  T#VMA>44ANri>l<*H  GO  TO  4 

Tf.C»r.f,T/v**l. 

noT«OOT*TfN 


*7- 


7 

IFIOrtI,u.ni*»|  r.o  to  7 

1PIRANOI  H>.GT,|nnT-ftT»!)/TCR  )  GO  T«  / 
•toTarwiT-ir.* 

«.«■  in  io  * 

FM>MANhW|S 1 

FviHAhiW  I  w  ) 

fcNO  ___________  __  __ 

*  'siiutirn  I  |nF  oomi  1  MAW,/  •«  ,tr } 
niMFNSlOA-  MAPI  1  ) 

LX.| 

HA«k 

1  «A.UAM(LX]*I  - 

H 

** 

FM.MiHN|V| 

F>04(||0l*FIMFV*FV»M*FM| 

IF  1 F |  ?.?,* 

IF  IF-I.J  «.<>.? 

FliBWWFl  AFlI/F 

Fv»V«Fl *FW/F 

CAII  THIMA.M4) 

ma*»ii.x|»um 

IF  ILX-Kl  7,7,4 

7  4F  Turn 

4  IX-IX-I 

HlaVA 

Fw«ur»i «Fu/f 

Ml  TO  •. 

III M A |  •0«S*|S.A4llA  —  Fll  | 

FNO 

VIMA)*O.S*< VA«VR-f V) 

Tui‘«Ol*nK-l-  TRS|M| 

P(MA)*0.4*(UA4UR-FW) 

COMMON  /VARS/  XO.YO.XM, YM.nw, VM,^M,  hip. 

|0(NA|>in{MA|«7 

COMMON  /STOPF/  NCI  ,  MrP,MAP|40),XI  *»0()(p  |.vt‘-000>,M(  5f»0i>  )  ,  W(  6i*o>i  t  , 

t<(MR)  «n.S*f|lA*MR«Fll) 

x  k/f  4000) .  ini  soon) 

V|MR)*n.S*(VA*VR*FV| 

XlM)»XM 

UJMRJ  *0.S*(WA*WfUFW1 

Y ( M ) ■ YM 

IM(MB) . iniNB)*? 

IMM)*IIH 

ITF«Y*I 

VI MJ  *VM 

IFIVRFI  .FO.O.)  WR|  TFIN.f.MFXI 

W|M ) *UM 

1  F  (  01*1  -OTMJ  4.  10,10 

IO|M>. IOM 

in 

n|rOT*Df>T*FKftA 

UPTURN 

MFIORN 

FNO 

FNO 

mWlliTIW  Imm.m 

FUNCTION  SOM 1 1.  « M,  X  ) 

COMMON  /STOmF/  NCL,FfP,MA**|  40  )  ,  X 1  4«00 1  ,  Y  |  *»000)  ,|l|  SOOO )  ,  V I  5n00|  , 

c 

SKIMMCR  |*j70 

x  wi  soon) .  in:  snnru 

OIMFNSION  XU) 

X(M|*X|N) 

a-xci  1 

Y|M)»Y(N) 

H-A 

IMM)all(N) 

on  4  |*|  ,M 

V|M)-VIM 

xm-xi n 

N|M)«M|N) 

IF  IXM-AJ  1,7,7 

|n(M|«|n|N) 

\ 

AaXM 

RF  TORN 

7 

FNO 

3 

A*XM 

WU*MMTRT-’frpinams) 

4 

CONTJNIIF 

C  SKIMMFR  1 V70 

SOM  * | R-AMIM-A) 

COMMON  /(.FOM/  PFI  ,01  (NS,OI  AMS,IHI  1 M«  M.0, SI  1 . XCftl • XCR7* XC R 3. kC A4, 

RFTIIRN 

X  SI  1 l,SLl7,XCHn,XCR|,XCRS,vrAl ,YCA?,YCH4,YCrt4,vrhS,VCFN,vrM7, 

FNO 

X  SlOl.Sl  07, SLOT, Si  04, SI  OS,  SLI4,S|I4,  MAn|N,01SN|,  XNJ.YNf. 

SORROOTINF  TRRXPfNARS) 

SK IMMFR  1*170 

COMMON  /T I MF /  T,nr,nTM,TFO,r>TSM»TFNntUTHP,TRFP 

COMMON  /STOKF/  Nf.l  ,MCP«MAP|30)«Xt 5000)  »  Yl 5000 )  *U|  5000  )  «  VI *000)  • 
X  W(SOOO). 10(9000) 


XNO,YN|l,  XN5.YNS,  V0II74) 

COMMON  /VARS/  XO,YO,XM,YM,OM, VM,uR, low 

COMMON  /fONS/  CH«SIA,MA*rM,CT»ACCnM.»'«e'>*CLMT(?MtTCnNM?S) 
COMMON  /HANO/  R,G 

COMMON  /HfRO/  nftPR,f)HPH,RPMF|  |40|  .MPkFUIO) 


1 

7 

3 

COMMON  /VARS/  Xn,YO,XM,VM,0M,VM,NM, |UM 

DATA  IC/1MR4/ 

M«  | 

MPCL*MAPINCt> 

IF  1 M-MPCI  »  4,4,3 

RFTIIRN 

3) 

37 

COMMON  /SW/  KSAM.KRFp 

U-OT 

CA1  1  irH<|SlOS,YCM6.XCHl,XCRS,>l  ,,r.S,IL) 

CAU  LCMM0,,YCP7,XCR5,XCRn,-l.,f  MU 

CAM  LCHKISLn,YCR3,0,,XCRO,U,F.N,TI  1 

CAM  if  MX  1  Si  I,YCR?,0.,XCHI,-1.,CI7.TLI 

4 

|0M* 101 M) 

C 

CMFCX  |  IP 

if  noM/tri  0,6,5 

IHUM)  1N.W.1 

5 

IO|M).|0M-|C 

\ 

IFIXO)  7.7  ,71 

7 

TTb-XO/IIM 

GO  TO  7 

in*lNv(tSIMAP,NCL«M) 

xn*x(M) 

VO»Y|M) 

IIMxIMM) 

-y**-yiMj _ 


WM-M(M) 

CALL  FRPINAPS) 
LM-LOCIXM.YMJ 

C10>CIWT(IM)/CLU1(|  0) 
ICLO-INTR(CLO) 


IHU-tJI  IN. 3, 3 
^  CALL  TWANMTT.O) 

1MVCB?-YN)  4,4, JH 
JFIVCM-YM) 

HP AM  FKOFIlF 


IFIU-TT)  |M.?7,77 
(.0  70|  )N,73),KSAM 
73  CAU  IR  ANSI  T  T  .0 ) 

CALL  SIFUO.,nRPK,RP«P|  ,50,1, 


C5*PU/-*»*- 

vw«<»«9vvii  *r  5<  vvr 

Wtwi 
VM  »-><■* 


PM.  ..... _ _  - 

<.n«minT  r  "  ifMK(<,«K.*i.Xll,*l,*.Tlt-M) 

rjiwr^M  /v# a*5/  xfi.Yp.yR.YM.oM.VM.Ni.,  in* 

|u|S  <i*PPftni |mc  ruff**;  MMFT«*-“  *  m»M  r(.MLf  f.nilinf^  WIT*-  iwt  5fr.MFN' 
fwliu  it  (it  km  OF  TKF  I  IMF  Y«5*X*»  .  II  15  A  PAMAMF T*-#  i'iM'.-AI|W.' 
INF  FWIM  V 


111  a  |l  *  5  1 
T*-/  Kali. 

TU«‘ 

Al 

(Ml  |  |  7.1  /«  1  7 


nnnmn 


C  |  ,  VM*V«**k  HJMM  -|5*HM)*#7 
C7»Yfi*VN-A|  *5*OM 
f.3*Yft*Yfl-AI  0  Al 

IFir.ll 

inr?>  /».  i  / « a» 


mi  in  i* 

e.?*r?/c  i 

r.w  3/r  \ 
nisfR*r  ?*f7-r.5 
|F|ni5fM|  >7. *7. <4 

Mf»»|5.5»MHlH5r,H| _ _ _ _ _ 

tmh.-w  ims-r.; 

IF  (1NAX.IT.0)  r.ll  in  17 
IF  (  TMM.I  T ,0)  Ml  1"  1  ? 

Tt«|M|M 
mi  to  i* 


Tf.JMAX 

IFIPl-TTI  17, H*l** 

1^  CAII  fMAMMIT,n» 

IFIXI.II.XM  ,AN|*.  XM.IT.XMI  1.11  TO  1#. 
IF  I  IT. I  I.TMAYJ  (.n  TO  17 
MFTOMN 


TI.FP  T*Of-T  T 
CAM  TMAM5ITT.I) 

MFTMHMI 

FNO _ _ _ 

•^TmmOMII  IMA  ANHK,SI 

5K|MMf-R  |*470 _ . _ 

- TUMMOM  /C0N5/  01.51*. 5|M.rM,Cl,A'‘r«*M.KMM.,riMTI7«i|,iri»»*M7M 

CfW*440M  /VAR5/  XM.Vn.XM.VM.IIM.VM.UM.lim 
MATA  lf/)A3K4/ 

COMMON  /KAMO/  H.r. 

WN,rT*5onit-Ai  ot.ihamoi  inm 

vt»r,M* » Awnc.ir.) _ - _ 

OMaC*VN-5*Vl 

VM,C*V1«5*VM 

WM»f.«*MAMMf.(r,| 

xo*xm 


5CAMfM  Fftw  | *‘T Fu  VAI  71  Il.7l.il  1 .1*1*1  I  WMFQF  Y  MCOlYV 

1*1 

J*S 

Vm{  l*.l)/7 
1FIX-/IM)  fa#," 


W*  Tn  <4 

laK 

IFt  1-1-1)  10.10,* 

F5‘ft  OF  U|MAttY  5FAMC*' 

7 1*JI  I) 


/  |  ,1*7.1-/  1 
X  I  •  C  X— 7  1 )  //  1 .1 
X.l*l/.l-X)/7l,t 
HaFII  | 

F  l-F(  i) 


l)F  la/  1  jaMF  I  ,1) 
f.|  allF  I  -?  ,*F  I 
(..JaMF  .1-7. *FJ 
Hl»3,*F l-OF | 
H.H3.*F.)-HF  J 


UFTIIWI 

FNII  _ . 

5t|5kOHl  l»»F  f/BO 

COMMON  /5T**RF/  NCI ,MCP, MAPI  50  |*XI  5001) )  ,V  1500ft)  ,Ol*fl00),V|  5000) 
X  WI5000I. 1015000) 

COMMON  /WAR 5/  yO-Vn.XM.YM.MH.  VM.OM.liih _ _ 


COMMON  /TIMF/  T.OT.ftTM.Tro.nTSM.TFMli.iiTMW.TMFM 
COMMON  /Ml  5T  /  FH|  10. 75I.FVI10, 75). I-WUO, 751. COM  10.7* 

JCFWl  10.  75  . Ml  75)  .  VVH I  75  )  ,WWft  I 

Ml  A5T»f 

on  *  i  »i.nr : 

L*I _ _ _ 


M|  A5T«WAPCI  ) 

15  (  /F»M|.  *5T*  *-MF  1*51 
IFII51/F)  5.4.7 
ItaMT/H  OATH  ^  I > 
nr  M  MaMF |P5T,M|  A5T 


VMaVI ►) 

WMaw|M) 

HTAr.aMnni  1ft|M).57)/7 
I  f-  (MfAM  1.1.7 

0 Al  I  5|F  SI-1. 5,0.3. M»|  ). I  ).1«.ftT.I»**l 

CALL  Sl£H--L^,0^,fVtl,»  - - 

r.AII  5IFT|-|.S,0.5,FM<  I  ,1  l.lO.ftT.WM) 

Ml  Tft  5 

IFIM1AI./7)  A. *. 5 

CINM.  ) afllin.il  |*m 

f  A|  |  5|F1|  — 1.5.0.  <,fMiU,l  »  .lO.ftT.MP) 

_ f  .Ml  5IF  Tl-l.s.0.5.f  FVI  I.IO.OT.WM _ _ 

r  M  I  5  |  F  1 1  —  I  .5,0«5,fFM|  |  ,1  ),|0.m.W*A) 
nf* ( i  }aitM|  i  )*liMan 
VM|I ) aVM| i ) ♦ vpan 
PM|l  )*MM|| 


VIlaYM 

ll)Mc  )  f  *  (  I  I'M/  If.) 

KFTIIHM 

FMM 

5IIRMOOT  IMF  HFFMXM.VM) 

COMMON  /VAW5/  XM.VM.XM.YH.IIMTVM.MM,  I  Dm _ 

MATA  If./ 143*4/ 

VMf|t|H*7.A  <  MMAXMf  VMt  YM| 

IIMallM-VMIIUMAXM 

VMaVM-VMORMaVM 

XU*XM 

_Ylla  Y** _ 

Iiimt  |rt|  lOM/iri 

MM  lllRM 

fMO  _  _ 

SMBWntlM^M  |MPM1 
5KIMMFU  1*4/0 
NlMOFIf 


COMMIV*  /f.fOM/  PF|,ni5M5,ft|AM5«lHl  IP.SL0.5l  |  «  XC  A 1 «  XCA7.  XT*  5.  <f  **&. 

X  51  1  1.51  17, Xr.Aft.Xf.nl, XCKS.YrMI  .YCM7.YCA3.VCMA.rr«5.VfM*.Yr«7. 

X  511*1.  SL07, 51.03,  51  OA.5I  05,  5113,51  K,  HA0|N,MI  5**| ,  XM).Y»*|. 

X  “NM.YMO.  XM5.VW5,  VO| |75| 
niMMOM  /VAM.5/  Xft.YP.XM  YM.MM,  VM,  HN,  I.IH 


COMMON  / MPRFI  /  7MPM< |0).FMPK(|0).OFMmk(1O) 
COMMON  / T  I MF /  T,ftT,ftTM. TFO.nTSM, TFNO.OTMP.TRFP 
COMMON  /mANo/  R.r. 

NIIMalMTMf  rilAOIt 
Oil  A  |«|.N!IM 


0051 1  ),ini*.||  ) ♦ilMvll*  *11 
ywMIi  >*wwm|i )«vM»VNvn 
WM  A  I  I  1  C  WH(<  I  ( 

H  fllMf  |N|IF 
|N||M 

- ‘itlDe- - - 

FM* 

\iinmiMl  |N|  M  1 1 1  i.f) 

COMPI1F  /VA** 5/  XO.Vn.XN.YP.ON.yM.w*,  It  ► 

rilMMfp-  /»!*-►/  T.OT.nlN.  l»n. <M55.  )*■*•**. 1-Jkm.  IMAP 

riiNPOM  /  5  f*  *m  F  /  Nfi  ,  nf  m*  .  M  A  P(  3n) ,  XI  5||0|«)  ,  V  1 5000)  *015000  t  ,  VI 5000  1 « 

X  Ml 5fH*f*i  .  1 1*1  Annul  _ 

Mf  aN API  Of I  )al 

) F  tU'-MfP)  3, 5, * 

1  WMlIf  I*.?) 

7  f  OWN  AT  I  I  fM  IKMAlif  nwi 

TFNt*.  T 
PM| 


3  f  A|  |  laiMN  |  NAP.Mf  .NTL.l  ) 
NaMAPlI ) 
fA|l  TM5IM) 

MF  1l»«M 

FM! 


HI  Af f OM'l  AT F  ^  •,|F|;i.  vAV  FO«  WA«|AM|F  w  HMFMF  X  J\  Tmi 

h|5TOi.wam  array  o»-  »»i ymmaioa-  n.  |m  a»i«-t)0  ,».Ti  w  .it.  A,m, 

/  )  5  Al«nn>  T«*  xial.  IF  4  F  Al  1 5  OnT5|i^  APR  AY  MASf.F  |T  fl*  TvJPOTF". 
Tl*  Tm*  NFAMF5T  f*"if  h|  |  . 
nlMm*!iii'  x|N> 


VO»M  Aft  |N*M  ANf>|  |«) 

I  FIT  A5p|  |  /NP*I  ,FWPM  .OFNPII  ,  1 0.  VP  ) -9  A*0L  I M  )  )  1,5,4 
HYPa5ilMTU*l5N|*nt5N|  »YIHYO) 

VX|,n|5N|/HYP 
VM I  a VO/HYP 

5IM5RAN 


IFI/-*.5*«  ANf*L  C  M  )  )  M.5.5 
I**/ 

VH«VM  1*51  3%if  ANflOlft  ) 

MP*  5|  naKAMnr.ir.) 

|  lIM.rt 


I  l.lainfi  I»>.YM) 

IF  1 1  1 .1.1  ¥  ,75 1  f.||  TO  4 

WMlfFlM.MIA) 

4 1 4  F|*MnAII  |m  .  *  HIRMM  |v  l«*f*) 
f.*l  1»*  »*|M 


IM  «  )  7,7.1 
IF  |N-h)  3.4.4 

5*1 

|.li  )n  4 
fih 

ln|<l* 7 _ 


M»*  1IIMM 


FUNCTION  TNTM(X) 
fllMWIf  /MAVlV  M 
1*1 Nf |X) 

IM  I  X-F*  OAT  II)  -CAAlfll  (Mil  7 


l»1*1 

TNJMal 

HFTliRN 

FA  O _ 

Fi'^’r  T  rn\  w  a'CVm.i  j  *  > 


YFL*lV 

YFf  *YFI  F-9 

IX.JV 
M AOUL  *VFI 
RFliiMN 


Function  manor ixs 
MANOR  •i,»RAHOI  lx  )-l. 
KFf'lRN 
FMl 

fUNf  TION  OMAx(X,l, 


Y*X<  l| 

on  i  io*( 

IF  I Y-x III)  7.1.1 

v»xm 

CONTINUE 


(■unction  i^nui 
iscn*o 

I F  (  a. 0F.0.  1  |Sf.N»| 
MfTUHH _ 


SORMOOT  |NF  nilk'» 

COMMON  /STORE/  Nf.l.  .**r.P,P  (  |  70950) 

COMMON  /RfRO/  0RP«,OMPH.RPRFL(50).HPkH.( 10)  , 

COMMON  'f.FOM/  PFI  *01  SN$,0I AMS*  TM|  I P.  SCO.  SI  I  •  XC* I  ♦  XCR7. *C«3.  xf  M4, 


X  <101  •  <1.02.  HO). si  04.SI05.  SH3.SU4,  HAMN.^I  SN| .  XN|.Y»I. 

X  XNO.VNO.  XN5.YN5.  VOL  I751.VMDI741  **  / 

COMMIT  /COWS/  Cli,SIA.SfR.f.R.CT.ACCOM.PRFS.CI/MT(?4).Tf«NS<?4) 

COMMON  /RTO|N/RR,Rf,,A,NST  / 

COMMON  /TIMF/  T.OT.UTM.  TFO.OTSM.  TFwn.l)IMP»T/*FP 
-  COMMON  /Ml5TZ  FO(10,Z5UFy.tl^25)-.FynOt75).CFOn0.2<^CFVnoT_73t._ 
If.FWI  10.7S).llR(75)«VR(751.Wftt75)«IHIR(?5l.VVM(?S)*MWR(74).rnw(.  (74) 
COMMON  /RANO/H.f,  / 

INIFCFM  MR.Hf,  ,R.f, 

namfi  isT/r.Fn/  pfi .oisws.nuMs.m  !P.sin.su.xcni.xci»7.xcM3.xrp4. 

X  SlINSI.I7.XC«n,XCRI.XCFS.YCHl.YCP7.,YCP<.YCR*.YCRS.YCMM.vCOr. 


x  si.oi.stfi7.si 0 -,stn4.st os.  sin. si  u,  rahim, nisi 

X  XNO.VNP.  XN5.YNS,  vm,  / 

NA*Ft  IST/CON/  CO.SlA.StH.CR.CT. ACCOM. MtffA.Ctuf 
NAMHIST/RFf,/  RR.Rf.,  A.NST.R.O  / 

NAMFI  IST/TIM/  I. OT.OIM, TFO.OTS**.  Tf  NO.I)Tj*P»TRFP 


MM ITFI A. | | 

_ I _ FORMAT ( *  I  *  .49X  .'SKI MMFM  1971'//* 

MR  I  IF  (fc.C.Fr  | 

MR? IF |4.mN| 

MR!tF(#>.RFO| 

PMllF(4.T|N) 

KT»f-IFO 

C  PREPARATION  OF  HI  MOOR  A  MS  FOR  OUTPUT 
THPR.O. 

on  7i  1*1. <0 

I  *10C1XCM-«I*0RPR  I 

IHI.I  1.11)  TRPfi.RPRFI  (  1  1/FIT*TRPR 


RPRFl  I  l|=RPRFU  I  J/(flMT(i)*FI  T*6. 7P3  JM*  (I-0.5)*0RPR) 
TMPR.O 

00  77  1*1.10 
I  «HlCIO,I*nnpRl 

iMPR.THPR.HPRmn/flMTIO/FI  T 

KPR.FL  ( 1 1  *HPRf  L.(  1 1  /(CLMTI  L  )*Fl  T*4.?R31«* f I-Q.5 1 *fi*PR I 
INI  *  1.  S7079*YCP7*YCR7*A*  ICONS!  Ip /VOl  (  M/Cl  MTU) 

0(1  7*  l«l.NCt 
ftmi>flt*C(.m'mu 
Oil  73  I*!. 10 
FO(  I.U*FII(  I.tl/FTMl 

FW( l,M«FH( |.( 1/fTwT 
fF'MF  ,1  )  -fFIK  I.ll/FIWT 
CFVII.I  )*CFV( I .LJ/FIMT 
fFWII.I  IrfFMM.l  l/FTWT 
COR Oil.  )*fOMMI  l/FTWT 


/R(U»VM(|  1/Fr.T 
WHI I  )«WR|I  | /FI T 
MOR(U«IMlR(t)/Fl  T 
VVRU  )  *VVR(  t  )/Fl  T 
74  MMRp )-WWR(l»/FLT 
_C  .PR  1NT  lN(«_n^_KI5TCK*«R  AMS  _ 

MKITH4.41 ) 

41  FORMAT  (  /  /«  30X  *  *  FLUX  PROFIIF  AT  OR  IF  ICfc  *  •//.  1  DISTANCE  FROM  AXIS*. 

X  4**X.  INMMRFR  PFR  UNIT  T1NM.//J 
f  A|  f  H|STI(HPRFL.l0«P.7.A0.TnTH« 

00  *7  J.1.10 

_S7 _ PIJ .  J)-.l»OMPR _ _ _ 

WMI  (*(4.41)  (  *P||,  II. 1*1.471. J*l  .101 

41  FORMAT! 101 F 10.4. 5X.40A) ,F 10.5/1 1 

MM  I  T  F  ( iS  «  44 1  T07M.THPR.TNI 

44  FORMAT!'  10TAI  •  .4RX.F  10.4//*  '  NOMRFR/llMF  THRU  OR  I F  |  f,F  *'.410.4, 

X  THEORETICAL  FIniRF  FOR  NO  SkJMmFk  INTERACTION  *'.F|P,4| 

_ MR  1 TF (4.441 _ 

44  FORMAT!  lM1.34X.»Rf  A*  FlUX  PROFILE  AT  EM  OF  F|  OmfTf  t'r# 

X  •  01  STANCE  FROM  AXIS' ,*0X, ‘NUMPFH  PFM  UNIT  TIMF*  .//l 
CAU  MISTI  (RPRFl  . 50. P.7.  100.  10TR1 
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XR'.flNF  AT  l  FAST  ONF  COLll«ION,  ANO  TMe  TMIRO  f.|VF5  IMF  SOM  n.  TmF  F 
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Xf M  HAH.'//) 

f A(  I  HIST3(P.FU.FV.FW.I  ) 
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RETURN _  I  . 

f  N II 

SHMMOUT |NF  RAM ( P. I • X 1 
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APPENDIX  3 

-  —7  •  >  1 

*  !  1 
No  Interaction  Beam  Intensity  t 

"  1 

We'!  derive  here  a  formula  for  the  beam  intensity  at  points  close  to 
the  skimmer. 

;  *  1 

!  '  1 

Let  D  be  che  n6zzle  skimmer  distance,  L  be  the  skimmer-observation 
point  distance  and  b  be  the  .skimmer  radius Consider  a  particle  in  the 
■skimmer  orifice  plane  at*  pqsition  (r  cos  0,  r  sin  0)  whose  velocity  has  a 
probability  density  function  proportioned  to  •,  , 

exp-[(u-l)2  +  (v-r  qos  9/D)2  +  (w-r  sin  0/D)2]/c2 

*  i  1 

That  is,  the  mean  velocity  vector  points  in  the  direction 

1 

■  (l»  r  cos  0,  r  sin  0)  '<  ■ 

1  D  D  : 


The  impact  position  of  the  particle  on  the  plane  normal  to  ,the  axis  distance 
L  downstream  of  the  skimmer  then  has  probability  distribution  1 

,  2'  ! 

— — -  1  [  (x-or1  cos- 0)2  +  (y-or  cos  0.)2] 


L2d2  e  Lc  ,  : 


where 


L  + 


.  If  the  flux  through  the  orifice  is  0  particles  per  unit  area  then 
the  flux  through  the  point  x,y  is  . 


-gTg"  e“(j~)2  [{x-dr  cbs  0)*  +  (y-dr  sin  0)^] 


0  (x,y)  =  0  Adr  /2lTd0 
J-i  •  o  o 

.  .  L  c  i 

This  integral  is  messy  in  general,  but  if  we  ask  only  for  the  centreline 
flux  (x  =  y  =  o)  we  obtain  •  ! 


,1  '  .b  ,  1  /I _ ,2  2" 2  . 

0r  .  =  /  rdr  J  d0  — 5-  -(r—J  u  r  1 
■  Lo  .  o  o  T  2  e,  Lc 
Lc  ir 


■which  one  can  easily  evaluate  to  obtain 


0Lo'“  (^W^L  +  D>]?>* 
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»moh 


